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Abstract: We determine the one-loop correction to the anisotropy factor for 
the square Symanzik improved lattice action, extracted from the finite vol- 
ume effective action for SU(AT) gauge theories in the background of a zero- 
momentum gauge field. The result is smaller by approximately a factor 3 
than the one-loop correction for the anisotropic Wilson action. We also com- 
ment on the Hamiltonian limit. 



1 Introduction 

Improved actions |2J have become a frequently used tool for doing Monte Carlo sim- 
ulations. Recently the square Symanzik improved action was introduced JJ, motivated 
by the desire to simplify perturbative calculations. From the numerical point of view this 
new improved action is not expected to be more optimal in removing lattice spacing errors, 
although simulations 0, suggest it is not performing much worse than the Luscher-Weisz 
choice of the improved action. 

The square Symanzik action was introduced to allow for a simple background covariant 
gauge condition. The background field calculation is particularly suited for computing the 
renormalized coupling constant, not only in the continuum ||, but also on the lattice with, 
or without anisotropy |7], ||. Improved anisotropic lattices are used both for thermody- 
namics || and for extracting glueball masses on very coarse lattices |§. In both cases the 
aim is to enhance the resolution in the time direction. Also for the square Symanzik action 
anisotropy was introduced and used in Monte Carlo simulations 0. This has motivated 
us to compute the one-loop correction to the anisotropy factor for this improved action, as 
it requires only a minor modification in the calculation already performed to compute the 
Lambda parameter for its isotropic version. 

We will perform the background field calculation for a finite volume at arbitrary 



anisotropy £, using the methods followed for the isotropic Wilson action |10| . Our results 
will include the Wilson action and the square Symanzik improved action. The one-loop 
correction to the anisotropy factor for the Wilson action was computed before by Karsch || 
in an infinite volume. Our finite volume calculation nicely confirms the universality of these 
results. 
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2 The anisotropic square Symanzik action 



For the Wilson action the anisotropy was introduced as follows ( *+p 



e SU(AO): 



yo x 



vc 1 E i^ + rtE^ 

i>j>0 iy^O 



P^ = 2KeTr(l-*Q 



x /I 



where (J3 = 2M / gl for SU(AT)) 



is required to guarantee that the symmetry of interchanging space and time is restored (in 
the infinite volume and continuum limit hence restoring Lorentz, or rather 0(4), invari- 
ance). As we have changed the discretization of the theory, also the Lambda parameter 
belonging to the running coupling will depend on the anisotropy parameter. Both 771 (£) 
and A(£) can be determined from a one-loop calculation. In this paper we will use the no- 
tation £(go) = V ^i^i 9o)£,i sometimes in the literature also denoted by 7. In the following 
£(y"o) wm be denoted by £ for short; from the context it should be clear when £ indicates 
the tree-level value. 

As was formulated in ref. ||, one can similarly introduce anisotropy for a tree- level 
improved action, 



S{{ci}) 



with 



E^TrE^{co(l-0+2 Cl 



C 4 1 



(3) 



^u = ^, & = £ 3 , & = £ 2 - (4) 

We wish to emphasize that the issue here is not to improve this action beyond tree-level. It 
would involve the extra non-planar Wilson loops that also appear in the isotropic case @. 
Its coefficients, as well as the one-loop corrections to c , C\ and c 4 will be doubled in number 
due to the anisotropy. After eliminating redundancies extra parameters will have to be 
determined, one of which can be related to r)\. However, additional physical quantities are 
required to fix all one-loop coefficients, making this a less than straightforward generaliza- 
tion from the isotropic case 
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The renormalization of the anisotropy parameter is, 
however, determined by requiring the restoration of the space-time symmetries in the con- 
tinuum limit, and can therefore be addressed without computing the one- loop corrections 
to the improvement coefficients. 

We impose the renormalization condition not directly by the requirement to restore 
the space-time symmetries, but rather by comparing the finite volume effective action in a 
zero-momentum background field derived from the anisotropic lattice action in eq. (3) with 
the result for the isotropic lattice action. We may also compare with the result obtained 
from dimensional regularization in the continuum. We will study the one-parameter family 
of actions defined by 

c = 1/(1 + 4z) 2 , ci = zc , c 4 = cl/c , (5) 

where z = corresponds to the Wilson action and z = —1/16 corresponds to the square 
Symanzik action, which is improved at tree-level to second order in the lattice spacing. 
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The relation C4 = c 2 /co allows for a simple background field covariant gauge condition 
that is easily generalized to the anisotropic case. 

?9f = >/^E^i i 1 + < 2 + A*)( 2 + 4*)) = o. (6) 

The background covariant derivative is given by D^(x) = f/ At (a;)$(x + /i)f/^(x) — and 
the quantum fluctuations are parametrized as U^(x) = e^^U^x) for a lattice background 
field Ufj,(x). The free propagators in this gauge (at U^x) = 1) are given by 

Ghost : P(k) = ' 



Ea6 (4sin 2 (A; A /2)+4zsin 2 k) 

Vector : /> (Jfe) = — , P ^^ v 0/7 , ^ . (7) 
* ^^(l + 4^cos 2 (V2)) 1 ^ 



3 Background field calculation 

We compute on a lattice of size iV 3 x 00 the effective action for a dynamical (i.e. time depen- 
dent) zero-momentum non-Abelian background field, Uj(x) =exp(c,-(t)) =exp(c"(t)T a /iV) 
and t)o(^) = 1 (the anti-hermitian generators T a are normalized as Tr (T a Tb) = —±5 a b). 
It is obtained by intergrating out all non-zero momentum modes. No integration over 
the zero-momentum quantum modes is included, which for a dynamical background field 
would lead to breakdown of the adiabatic approximation near c = 0, where the classical 
potential is quartic |T3(]. We will follow closely the methods developed for the isotropic 



Wilson action, described at great length before |JT0 |. The effective action is given by 

where a t = L/^N is the lattice spacing in the time direction, L the physical size of the 
volume, Ft = e a b e c^ the field strength and Vi(c) is by definition the rest of the effective 
potential. All that is relevant to know is that at 0(c 4 ) its coefficients are fixed uniquely 
by an abelian background field, unambiguously separating (F^) 2 from Vi(c). We have 
ignored terms that vanish in the continuum limit. Furthermore, the renormalization group 
to one-loop order implies g^ 2 = — llA/'log(a s A)/247r 2 , where a s = L/N = a 4 £ is the lattice 
spacing in the space directions. We have therefore introduced the renormalized coupling 
g- 2 = g 2 -llAf\og(N)/2ATT 2 = -llA/'log(LA)/247r 2 . 

We note that A, 771, a\ and a.2 depend on £ and z. Universality requires that physical 
quantities, as well as the background field effective action, are independent of these pa- 
rameters in the continuum limit. This implies that both (a\ — llA/'log(LA)/247r 2 — 771 /27V) 
and (a 2 — HA/'log(LA)/247r 2 -|-?7 1 /2A/') are independent of £ and z. For isotropic lattices 
(£ = 1), for which r]i = 0, this implies that 

a 1 -a c 1 = a 2 -a c 2 = llAf\og(A/A c )/247r 2 (£ = 1), (9) 

where a\ are the values for a fixed isotropic regularization, like dimensional regularization in 
the continuum, or the isotropic Wilson action. For anisotropic lattices we follow Karsch || 
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by defining 

c T (0 = ^(1) c CT (0 = a 2 (l) - « 2 (0 (10) 

and one easily derives 

77! (0 = -A/-MO - c T (0), A(0M(1) = exp (-127r 2 [ C(T (0 + c T (0]/l W) . (11) 
We can summarize these various relations between the Lambda parameters also as 

A(£, z)/A(?, z>) = exp (l2n 2 [ ai fo z) + a 2 & z) - a^', z') - a 2 (?, z')]/llM) . (12) 

4 Analytic results 

The coefficients ot\ and a 2 are determined by working out the determinants of the quadratic 
fluctuation operator. This section can be skipped when one is interested in the numerical 
results only. To all orders in the background field and to quadratic order in the quantum 
field, using eg = and c^x + p) = cfa) = c£(f), we find (k, i e {1, 2}) 

S * = ^E Tr ( E ^ kU ' 2 {(Si^{t)-2) {D k ^{x)-D io q k ^x)f - (13) 
^W^( t ){[^( x ) + ^^ fc A( ;r )^fcA( a; ) + ^A^(^)] + [qa>(x),q k ji(x)] + 

(l+k-2)D e >[Ux), DM*)]}} ~ 4 {E ^UWf) • 

We have introduced the following convenient shorthand notations 

Skfi/uit) = 2 — Skjx/u(t) — Sl^ ec ,(t), S kjiii> {t) = Sl^ ei/ (t) — Skp,,eu{t), 
SkMu(t) = e k ^e ee » {t+k ®e- kd ^ t+eo) e- e£ » {t \ (14) 

as well as "doubled" covariant derivatives and quantum fields (qifi(x) = q ll (x) and Di^D^) 

D 2 ^(x) = e 2£ ^$(x + 2fi)e- 2 ^ - q 2fl (x) = 2q^(x) + D^(x). (15) 

Therefore D 2jl = 2_D M + D 2 and the gauge fixing functional can be written as 

= V^12^ zk ~ lD Um(x). (i6) 

S 2 was obtained by adding — X^Tr fig{x)/ g 1 to the action S. Under an infinitesimal gauge 
transformation, q k jx{x) transforms as D k pQ{x) for k = 1,2, giving for the ghost operator 

M^v^E^^Aa- (17) 

H,k 

It is now straightforward to compute the functional determinants. Simplifications can 
be made due to the fact that we can split the one-loop correction in a purely kinetic part 
for which we can drop all terms of higher than second order in c and a potential term for 
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which the time dependence of the background field can be ignored. We find the following 
results 



a s+1 (X) = a^ 2) {X) + ?aW>{X) + + a^(X) 7 



S s > 2 h 



with X = {£, z, TV, N}, s = 0, 1 and 



a 



a 



a 



(*) 



117V, , An AT 
2^ l0g(iV) + 6X^ 

AT 



E 



8vr^ 3 iV 3 



E 



12tt£/A 3 



fc^O 



dfc |2(4-3 Pl/ )( 
<9 2 



3f 



+ 



<9 



4 



6CP 



dk (d^d v 

Af 2 
Cp" + 4(1 



(19) 



, 2 d 2 logP 



+ 4P 



1 + 4^, P = P(fc, £, z) is the (rescaled) propagator, and d^, p^ and g M are simple 



where £ 

trigonometric functions (momenta are given as fc = (/co, k) 

-l 



4 sin 



t/eo) (l+4z 



cos 



k ))+uj x 



(k ,2irn/N), m E Z N 
r 2 E 4sin2 (5^) (l+4^cos 2 (iA; 



dp = cos(/c M ) +4z cos(2/c /j 



w , p At =[l+22;(l+cos(A; M ))] ; g M =jy 

We note that in the continuum limit, N — > oo, the total derivative terms will only 
get contributions from near k = 0. It can be shown that 



-i 



;i+cos(/g)(2-c?v) 2 .(2o) 



lim a ( b s ' 2) {X) =M 



N^oo 



S + l 

72tt 2 



+ s 



804 



457T 2 



a rl 



-(4tt)" 2 - 0.619331710- •• (21) 



where a± is a constant introduced in ref. In particular these boundary contributions 
are independent of £ and 2 and drop out in the computation of the quantities in eqs. (10- 
12). The remaining terms can be converted to finite integrals in the continuum limit, 
replacing iV~ 3 by (2tt)- 3 d 3 k. 

For isotropic (£ = 1) actions, a^ 2 '^ = a^' 2 ^ and a^' 2 ^ are independent of s. Conse- 
quently aii — a 2 = A/"(l/247r 2 — 8a 4 )/5 is independent of the regularization employed, cmp. 
eq. (9). The non- vanishing value of this difference is a manifestation of the breakdown of 
Lorentz invariance in a finite physical volume. At finite N our results in eqs. (18-19) are 
exact. All integrals over ko can be performed analytically, yielding sums over the N 3 - 1 
non-zero spatial momenta of explicit analytic expressions in z, £, M and k, which are read- 
ily evaluated numerically. It can be shown that as an expansion in 1/N, terms linear and 
for z = —1/16 (i.e. with improvement) quadratic in the lattice spacing are absent. For 
z — 0, where p M = £ = 1, q^ — (1 + d^) and d^ = cos(fc M ), dramatic simplifications occur. 
There are in particular for z = more efficient ways to compute the coefficients, but we 
will not dwell on this any further. 

We list the following k integrals required to evaluate the expressions in eq. (19) (y(u) 
is defined by Cj = 2 sinh(±?/) and m > 0) 

m— 1 ( 



27T J — 7T 



dfco P" 



<9 



1 



(m-1)! 



del; 2 



2sinh(?/)vl + 4l(F 



5 



5 1 



1 + Aze-y + VI + 4ze-y 



^ = ^L dkoP ° p =^Ty\-d^) i^^i + — ^ — /' 

5 Numerical results 

By extrapolating in 1/N the explicit expressions for aci and «2 on a finite lattice, we can 
extract their continuum limit to at least nine digit accuracy (we evaluate the momentum 
sums on lattices with N = 3 to N = 99). Our results are valid for arbitrary SU(A/"), 
where 011,2 can be written as aAf + b/Af. Using eqs. (10-11) we reproduce results obtained 
by Karsch || for the Wilson action. The ratio of the square Symanzik action Lambda 
parameter to the Wilson action Lambda parameter is obtained using either eq. (9) or 
eq. (12) with (2, £) = (—1/16,1) and (2', £') = (0,1). The result was already reported in 
ref. 0. One finds 

/ 4.0919901(1) for AT = 2, 
sq/ w ~ \ 5.2089503(1) for Af = 3, 1 ' 

agreeing with two alternative recent determinations based on the heavy-quark potential 
and twisted finite volume spectroscopy [fT^| . 

In table 1 we give results for the square Symanzik action at some selected values of £ 
between 1 and 20, likely to be of use in simulations, as well as for the Hamiltonian limit, 
£ = 00. The £ dependence is illustrated in figure 1. Results for the Wilson action are given 
for comparison. We only present T]i and A(£)/A(l) for SU(2) and for SU(3), since one can 
use eq. (11) to extract the values of c T and c a . These can furthermore be used to extract 
the results for any other number of colors, since 

CtAM) = 22^,(2) + ?^W(3). (24) 

We note that in all cases the value of T]i is reduced by a factor of approximately 3 for 
the improved square Symanzik action as compared to the result for the Wilson action. 
Indeed for the simulations performed in ref. [H] a reduction with approximately a factor 
2.5 for the measured value of 771 can be deduced (whereas the results obtained from the 
Liischer- Weisz and square Symanzik action agree within errors) . We extracted rjx using the 
one-loop truncation of eq. (2). At the rather strong coupling employed in these simulations 
the measured values of r\\ themselves should of course not be expected to agree with the 



perturbative results |L4| 
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6 Hamiltonian limit 



Here we briefly discuss an interesting feature of the Hamiltonian limit, i.e. £ — ► oo. It turns 
out that the potential Vi(c) has field dependent contributions that diverge in this limit. 
This divergence, however, only occurs for z 7^ 0, in particular for the improved square 
Symanzik action, i.e. at z — —1/16. At first sight this may seem puzzling. However one 
should notice that when modifying the action also the measure of integration has to be 
corrected. Such a modification of the measure can of course be absorbed in the action, as 
is usually done and is of the same order as the one-loop corrections, giving rise to a term 
J2a 5V{. c i{t))- Since the effective potential appears in the action as J2t a tVi{c{t)), with a t = 
L/N£, we conclude that the total contribution to the effective potential, due to correcting 
for the measure, is linear in £ (vanishing for N — > 00) and given by iV£ J2i SV(ci)/L. It is 
hence much more natural to redirect any terms linear in £ to the measure. 

To determine 5V we compute Vi(c) by taking an abelian background field, suitably 
extended to the non-Abelian sector. For SU(2) this extension is achieved by substituting 
Ci = VE« c t c i m the Abelian background link variable Uj = exp(±iCj<r 3 /N). This can 
be generalized to arbitrary gauge groups following the methods described in ref. fll5| , but 
we will for the sake of presentation only consider the effective potential for SU(2). At 
0(c 6 ) there are additional terms that vanish for Abelian background fields, but they do 
not concern us here (and have finite limits as £ —>■ 00). Along the lines described in ref. || 
anisotropy is easily incorporated and one finds Vi(c) = Vf°(C) — V^ ab (0), where 




1 



4UI 



u 2 UJ 




(25) 

with Xjin^Cj) = 1 +4,2 cos 2 ((vm,- + \Cj)/N) and u 2 {n,C) = 4Aj sin 2 ((vm j + \Cj)/N) 
(a more detailed derivation for the isotropic case will appear elsewhere |[16|| ). One can 
verify that this gives the correct continuum limit at fixed £. For the Wilson action (.2 = 0) 
one finds V£ b (G) = L -1 iV£X^ /54asinh(a;/2£) (the apparent divergence for z = can be 
shown to be field independent). At finite N we find 5V(ci) = J2n^o^°S {^i( n ii Ci)/\( n ii 0)). 
The remainder, denoted by Vi(c), can easily be shown to have a finite limit for £ — >• 00. 
We conclude that the Haar-measure, dUi(c) = N~ 1 (27rC i y 2 sin 2 (Ci/2N) H a dc°;, is to be 
corrected with a factor exp [— 5V(ci)], giving the exact measure at finite iV for deriving 
the effective Hamiltonian from improved actions. It is not too hard to show that up 
to exponential corrections in N, we have exp[— 5V(ci)] = [1 + 4z cos 2 (Cj/2iV)] /(l + 4z). 
Indeed at z = —1/16 the rescaled measure is flat to 0(c 3 /iV 3 ). 

It is also interesting to point out that one finds LVi(c) = 7i(c") 2 + C(c 4 ), with = 
7^ — 2z/N^yJl + Az + 0(N~ 3 ). Since at z = —1/16 on-shell improvement should imply 
that spectral quantities have no 0(1/N 2 ) errors, the 1/N correction to 71 can be removed 



by a non-local field redefinition, as was explained to some detail in ref. [T^j. The field 
redefinition is designed to remove the next-to-nearest couplings in the time direction (not 
listed in eq. (8) since they are irrelevant in the continuum limit). This non-local effect 
disappears in the Hamiltonian limit, as it should. 
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7 Conclusions 



We have calculated the one-loop correction to the anisotropy parameter for the square 
Symanzik action, using a zero-momentum background field calculation in a finite periodic 
volume. For the Wilson case we retrieved the results in an infinite volume of Karsch ||, 
which is a rather non-trivial check of universality, since even in the continuum the periodic 
boundary conditions break the Lorentz invariance. We find that the size of the one-loop 
correction to the anisotropy is reduced, both for SU(2) and SU(3), by approximately a 
factor 3 when using the square Symanzik improved action instead of the Wilson action. 

We have also discussed the "Hamiltonian limit" of the zero-momentum effective theory, 
where the lattice spacing in the time direction is reduced to zero. We show how the 
integration measure is to be improved, defining the inner product on the Hilbert space 
involved in extracting an effective Hamiltonian from the effective action 
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Table 1: Results for rji, the one-loop correction to the anisotropy £, and the 
Lambda ratios for SU(2) and SU(3) Wilson and square Symanzik improved 
lattice actions. 
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Figure 1: Comparison between the square Symanzik (sq) and Wilson (W) action results 
for and A(f)/A(l). Bottom figure for SU(2) and top for SU(3). 
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